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ABSTRACT 

We report the first localization of diffuse, non-thermal, X-ray emission from a nearby galaxy cluster. 
Using Chandra data, we have isolated a diffuse non-thermal X-ray component with a photon index, 
Y x = 2.21±o;*5 and a flux of 9.5^1^ X 10~ 5 photons cm" 2 s" 1 keV" 1 at 1 keV, that extends from 
~1'.5 to ~2'. 5 to the south of the X-ray flux peak. Comparison to simulations implies that the diffuse 
non-thermal emission is produced by primary electrons, accelerated at shocks to relativistic velocities. 
Using these results and the flux and hardness maps produced with data from the Chandra Advanced 
CCD Imaging Spectrometer, we conclude that a smaller subclump coming from the north merged 
with ICI262. The offset of the cD galaxy from the X-ray peak and large peculiar velocity indicate 
that the subclump's impact parameter was to the west and on the near side of IC1262. 
Subject headings: galaxies: clusters: individual (ICI262) - galaxies: intergalactic medium - galaxies: 
X-rays - galaxies: non-thermal emission 



1. INTRODUCTION 

Recent Chandra observations indicate that rela- 
tivistic plasma injected into the intracluster medium 
(ICM) from radio sources eventually detaches from 
the radio so urce, forming bubble s of radio p lasma 
in the ICM llFabian et al. I 12003 iHeinz et all 120021: 

UD- I 



McNamar a et al. 1 1200 111 . In the model proposed by 



Enfilin (1999), these radio ghosts survive in the ICM, 
and provide a seed population of Cosmic Rays(CRs). 
Diffuse non-thermal (NT) emission is produced when 
merger induced shocks re-accelerate, via the first order 
Fermi process, this seed CR population. Current evi- 
dence suggests that these radio ghosts contain a signifi- 
cant population of proton s as well as electrons |EnBlin 
1999). Since lBlasili|2001ft demonstrated that diffuse NT 
X-ray emission could be produced by either primary elec- 
trons directly accelerated at shock fronts, or secondary 
electrons produced during proton-proton collisions, there 
are two possible sources for the observed diffuse NT emis- 
sion. To determine the mechanism that produces diffuse 
NT emission requires accurate measurement of the spec- 
trum and locatio n of the NT emission . 

Simulations bv lMiniati et al. 1 l)2001j) show that diffuse 
NT emission occurs in a broad mass range of clusters 
with a luminosity proportional to the X-ray temperature, 
making cool clusters and groups an important diagnostic 
for understanding which population of electrons produces 
diffuse NT emission. They find that spectral index of the 
diffuse NT emission is dependent on the electron popula- 
tion producing the emission, such that the spectral index 
of diffuse NT emission produced from primary electrons 
has a steep spectral index (aj c >1 .1), while for secondary 



it is systematically flatter (a ir <1.1) 
Hud son. Henriksen. fc Colafrancesco I l)2003fl reported 

1 Laboratory for High energy Astrophysics, NASA/GSFC 



detection of diffuse NT X-ray and radio emission from 
ICI262, a poor cluster of galaxies. The X-ray detection 
was made using the BeppoSAX Medium Energy Con- 
centrator Spectrometer (MECS) detector, and the radio 
using the NRAO VLA Sky Survey (NVSS) and West- 
erbork Northern Sky Survey (WENSS). Although the 
MECS was able to constrain the region of the NT emis- 
sion, it does not have the spatial resolution to identify 
the source of the NT emission. In this paper, we show 
that the Chandra Advanced CCD Imaging Spectrometer 
(ACIS) has the energy range and spatial capability to de- 
tect, localize, and characterize the spectrum of NT X-ray 
emission from low temperature clusters. These results 
are compared to simulations of cosmic-ray acceleration 
at shock fronts. 

Throughout this letter, we assume a Hubble Constant 
of H = 65 /i65 km s _1 Mpc -1 and q = h. Quoted 
confidence intervals are at a 90% level, unless otherwise 
specified. 

2. OBSERVATIONS AND METHODS 

The IC1262 galaxy cluster is a poor cluster of 
galax ies located at (J2000 ) 17 h 33 m 01.0 s , +43° 45' 
28" labeling et al. I 1199^1 with a redshift of 0.0343 
l)Colless et al. Il999|) . so that V = 46 hgg 1 kpc. It was ob- 
served by the Chandra ACIS S3 chip on 23 August 2001 
for ~ 31 ksec. The total count rate of the uncleaned 
data is 10.2 ± 0.02 cts s _1 with a peak of 37 ± 3 cts 
s _1 , which is well below the saturation limit of the ACIS 
chips operating in Faint Mode (170 cts s _1 ). Figur&IDis 
a color coded intensity map that shows the full 8'.5x 8'. 5 
image in the 0.3-8.0 keV band. The image was obtained 
using processing routines outlined in the CIAO 2.3 Sci- 
ence Threads 2 . The data was CTI corrected and cleaned 
for flares, point sources, and anomalous high background. 

2 http : //exc . harvard . edu/ciao/threads/ index . html 



2 



Hudson, and Henriksen 




Fig. 1.— Chandra ACTS Flux Map of IC1262. The white figures Fig. 3.— Logscale Hardness Map of IC1262. (2.0-8.0 keV)/(0.3- 

are regions extracted to search for possible non-thermal emission. o OkeV") 
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Fig. 2. — Temperature Map of IC1262. Regions with a (*) are 
have Temperatures given for the double thermal model. 



Exposure maps were created for 17 different energy bands 
to ensure proper exposure correction. Background was 
taken from the CALDB 2.21 blank-sky datasets. 

To construct the temperature map (see Figure - } we 
divided the chip into quadrants, north, south, east, and 
west. Spectra were extracted from pie shaped regions 
with radii differences such that there were at least 20000 
Counts in each region. In order to ensure proper back- 
ground subtraction, we eliminated periods of high back- 
ground, following the Filtering Lightcurves thread 3 . In 
order to account for background variability, the back- 
groun d was normalized to the source in the 10-12 keV 
range l)Markevtich et ai~ll200.lt) . The data were grouped 
so that there were at least 30 counts per channel. All 
models included either an acisabs component or had an 
acisabs corrected ancillary response file in order to ac- 
count for the time dependent absorption (at low ener- 
gies) of the ACIS window. Results are given in Table - ] 

3. RESULTS AND ANALYSIS 

From the flux and hardness map, Figure - ]and Figure- 
~ ] respectively, we chose four regions that have a strong 
gradient in flux and hardness. These regions are shown 
in Figure - ] We considered four models of background 
subtraction, to ensure any NT detection would not be 

3 http : //cxc . harvard . edu/ciao/threads/f ilter_ltcrv/ 



the result of undersubtracted particle background. These 
four background subtraction techniques were: (1) stan- 
dard blank-sky fields extracted from the same region on 
the S3 Chip as the source a nd normalizing in t he 10 - 
12 keV range as described in iMarkevtich et, al. I <|2003|) . 

(2) extracting a region far to the north of the peak flux, 
modeling it and then using the model times a free con- 
stant in our source + background spectrum modeling, 

(3) freezing the background model described in (2) with 
a constant proportional to the difference in size of the 
source and background region, and (4) subtracting the 
size normalized background spectrum from the source 
spectrum. In all cases the source model parameters are 
consistent within 90% confidence, indicating that all four 
methods produce the same results. 

We modeled our four "possible shock" regions with 
a single thermal component model, a double compo- 
nent thermal model, and a thermal model with pow- 
erlaw. The addition of the powerlaw component does 
not significantly improve the fit in the East and West 
Regions. In addition, Figure - ] shows that a thermal 
model fits the entire spectrum in the East Region. In the 
case of the North Region, adding either a second ther- 
mal component, or a power-law component improves the 
fit. The single thermal model fit to the South Region 
has clear high end residuals (see Figure - }. Although 
the addition of a second second thermal component im- 
proves the fit (Ax 2 ~ 16 for a reduction of 2 degrees 
of freedom(DOF)), the addition of a powerlaw compo- 
nent provides the best fit (A% 2 ~ 23 for a reduction of 2 
DOF) and reduces the residuals (see Figure - ]) . Although 
we cannot rule out a thermal interpretation of the sec- 
ond component, we argue a non-thermal interpretation 
is better statistically and physically. The addition of a 
powerlaw (37% probability of exceeding % 2 = 74 for 71 
DOF) produces a better fit than adding a second ther- 
mal component (18% probability of exceeding x 2 = 81.8 
for the 71 DOF). The 90% lower confidence of ~7 keV 
implies a shock with a Mach nu mber of ~3, which is 
higher than observed in clusters l)Cabici & Blasi ifeOOSl 
and references therein). 

In order to confirm that the nonthermal detection was 
not simply undersubtracted particle background, we re- 
peated our four different techniques for handling the 
background that we had done for a single thermal model. 



Cosmic-ray Acceleration at the Shock front in IC1262 
Table 1. Spectral Fits of Regions. 
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Region 


Modelt 


kTi 


Ni 


Apec 


1.27 


Ni 


Apec+Apec 


1.88 


N 2 


Apec 


1.34 




Apec 


1.57 


N 4 


Apec 


1.88 


N 5 


Apec 


2.09 


Wi 


Apec 


1.80 


Wi 


Apec+Apec 


1.98 


w 2 


Apec 


2.09 


w 2 


Apec+Apec 


2.38 


w 3 


Apec 


1.93 


Si 


Apec 


1.62 


Si 


Apec+Apec 


2.54 


s 2 


Apec 


1.60 


Ei 


Apec 


1.62 


Ei 


Apec+Apec 


2.08 


E2 


Apec 


2.02 


E L » 


Apec+Apec 


2.30 


E 3 


Apec 


2.37 



Abundance kT 2 (keV) 



X 2 /dof 



+0.03 
-0.03 
+0.17 
-0.15 
+0.12 
-0.04 
+0.13 
-0.16 
+0.43 
-0.37 
+0.95 
-0.58 

+0.09 
-0.12 
+0.12 
-0.13 
+0.18 
-0.12 
+0.22 
-0.21 
+0.16 
-0.16 

+0.05 
-0.05 
+0.26 
-0.22 
+ 0.08 
-0.10 

+0.05 
-0.06 
+0.32 
-0.16 
+0.21 
-0.18 
+0.29 
-0.36 
+0.37 
-0.37 



0.14^ 



,+0.03 
-0.03 
44+ ' 17 
,+0.05 



0.14 
0.12 
0.13 
0.12 



0.03 
+0.07 
0.05 
+0.14 
0.08 
+0.34 
0.12 



0.25 1 



-+0.06 

-0.05 

0-38_ ; 09 

33+ 11 
u,oo -0.08 
n c 7 +0.25 
u -°'-0.16 

U - 1S -0.06 



0.27"* 



7+0.05 
-0.04 

0.63+ 013 



0.17 

0.22 
0.49 
0.23 
0.40 
0.39 



0.22 

+0.06 
0.05 

+0.05 
0.04 
0.28 
0.13 
+0.11 
0.08 
+0.23 
0.20 
+0.25 
0.17 



0.83 



+0.16 
-0.04 



0.61 
0.28 

1.05 

0.82 
0.23, 



+0.20 



+0.09 



+ 0.21 
0.07 



+ 0.21 
0.08 



+0.09 
0.21 



173.0422/98 
134.0729/96 
119.7284/91 
57.09452/76 
65.19046/59 
55.56707/48 

126.3223/108 
113.5917/106 
136.9889/101 
123.7163/99 
118.0369/107 

203.8204/103 
147.7749/101 
134.1819/97 

114.6884/101 
87.91611/99 
90.03336/89 
83.43614/87 
109.8093/99 



f All models include photoelectric absorption using Wisconsin cross-sections frozen at galactic value of 2.46 X 10 20 cm 2 . 



As in the case with the single thermal model, the results 
are consistent within 90% for all the background tech- 
niques used (see Table-0). Also, we point out that the 
lack of a NT detection in the West and East Regions pro- 
vides more evidence that the detection in the North and 
South Regions is not simply particle background, since 
that would be observed isotropically across the detector. 

The radio halo detect ed by 

iHudson. Henriksen. fc Colafrancescol l)2003|) appears to 
extend to the south along the elongated X-ray isophotes 
and through the South region (see Figure-OJ). However 
the ~1'.0 resolution of the NVSS and WENSS precludes 
the removal of radio point sources. We, therefore, con- 
sider two possibilities: (1) that the extended southern 
radio emission is dominated by radio point sources, 
and (2) that it is diffuse radio emission. In the case 
that the radio emission is dominated by radio point 
sources, we argue that this would have no affect on our 
X-ray results. Relativistic electrons diffusing from a 
point source at their Alfven speed, would only reach 
~ 5 h^" 5 kpc (~ 6" ) (assuming a 1 /iG magnetic field) 
in their radiative lifetime (~I0 8 years). Therefore any 
radio point source with detectable X-ray emission would 
appear point-like and would have been removed from the 
Chandra ACIS data. In the latter case, the diffuse radio 
emission is emitted by the same electrons producing the 
inverse compton emission ( as in the invers e compton 
scattering-syntrotron model ijRephaeli lll9 79 , ,. We also 
note in passing that there are several identified radio 
point sources (including the cD galaxy) in IC1262, which 
could provide the seed population of CRs necessary in 
the model proposed bv lEnBlinl l)1999[) . 

4. DISCUSSION 
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Fig. 4. — Spectral fits to regions of compact isophotes: top is 
single thermal fit to East Region, middle is single thermal fit to 
South Region, and bottom is thermal + powerlaw fit to the South 
Region. 



A comparison of Figure-^ and Figure-|2l to flux and 
temperature ma ps created from simulations of c luster 
mergers done bv lRoettiger. Loken. fc Burns I l)1997|, and 
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Table 2. Nonthermal detection in the North and South regions. 



Region Background Model kTi (keV) Abundance 


Tx 


Non-thermal 


X 2 /dc.f 






Normalization(x 10 — 4 ) 





North 



blank-sky 
model-bkg 
model-bkg(Frzn) 
-local region 



1.32 
1.29 
1.31 
1.30 



+0.15 
0.12 
+0.16 
0.13 
+0.17 
0.14 
+0.17 
0.14 



> l-00±o >0 79° 

> 1.00t >?3° 
0.39i>°f 
0.43±>°f 



2.36 
2.30 
2.41 
2.40 



+0.14 
0.12 
+0.19 
0.17 
+0.24 
0.45 
+0.23 
0.42 



0.79 
0.78 
0.56 
0.58 



+0.09 
0.34 
+0.11 
0.35 
+0.22 
0.35 
+0.21 
0.36 



52.93678/50 
54.05200/49 
57.22508/50 
54.61599/50 



South 



blank-sky 
model-bkg 
model-bkg(Frzn) 
-local region 



1.44 
1.40 
1.44 
1.43 



+0.18 
0.13 
+0.18 
0.14 
+0.20 
0.15 
+ 0.19 
0.14 



> 1 00 +>a0 

> i-uu_ 62 
n 4q+>0-51 
u - 4y -0.32 

> 1 oo +>0 -° 



2.21 
2.14 
2.33 
2.27 



+0.14 
0.15 
+0.16 
0.17 
+0.33 
0.51 
+0.18 
0.38 



0.95 
0.91 
0.59 
0.78 



+0.11 

0.25 
+ 0.13 
0.27 
+0.26 
0.42 
+0.11 
0.46 



74.33978/71 
77.66279/70 
90.72980/71 
76.28895/71 
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Fig. 5. — Digital Sky Survey Image of IC1262 with Chandra 
ACIS X-ray flux contours (blue) and NVSS radio flux contours 
(red). The X-ray contours are 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.75, 
1, and 1.3 X 10~ 3 photons cm -2 s — 1 arcmin - 2 . Radio contours 
are 0.45 1.35, 2.70, 4.50, 6.30, 9.00, 13.50, and 15.75 mjy beam" 1 . 
The black and white crosses are Faint Radio Sources at Twenty 
centimeters. 



iTakizawal l)2000|) . indicates a recent (~0.25 - 1 Gyr 
ago) north-south merger of a smaller subclump (mass 
ratio ~1:2 - 1:4) with IC1262. These simulations show 
the general characteristics of the temperature and mor- 
phological structure: elongation of the isophotes to the 
south, compression to the north, and the hottest gas 
south of the center. Since cD galaxies in relaxed sys- 
tems will reside at the bottom of the potential well (e.g. 
lOeeerle fc Hill 1 1)200 1)) ). the offset of cD galaxy from the 
X-ray peak and its high peculiar velocity (Av = 453 km 
s _1 ) provide additional evidence of a merger. The kick 
to initiate sloshing is paral lel to the impact parameter 
ijTittlev fc Henriksen 1 12003) . As the simulations indi- 
cate a recent collision, the subclump must have passed 
to the west, pulling the cD galaxy to the west. 

The elongation of the X-ray isophotes toward the south 
that overlap the diffuse radio emission open up the possi- 
bility that the elongation is due to non-thermal emission 



rather than a merger. However, the evidence for a merger 
comes not only from the similarity of the southern X-ray 
isophotes to those seen in the simulations cited above, 
but also from the compressed northern X-ray isophotes 
and the very hot region that is south of center. 

Comparing our re sults with the simulations done by 
IMiniati et al. 1 1)2001)) . we argue that the localized diffuse 
NT emission is from a population of primary electrons. 
IMiniati et al. I l)2001l) find that primary electron reaccel- 
eration follows the shock region closely. In the central 
region (r < 0.3 Mpc), they find the greatest difference 
between the spectral index produced from primary elec- 
trons (ai C >1.1) and secondary electrons (aj c <1-1). Our 
90% confidence photon index Tx = 2.07 - 2.36 which cor- 
responds to a spectral index a^ c = 1.07 - 1.36 is consis- 
tent with the spectral index that IMiniati et al. I 1)20011) 
found for primary ele ctrons. Using the best fit Fi c -T x of 
IMiniati et al~l (f2001^ . we find a flux of (10 - 48) x 10~ 13 
ergs cm~ 2 s _1 over 0.13 to 100 keV for primary electrons 
in a 1.94 keV cluster (IC1262's temperature as dete r- 
mined by Hudson, He nriksen. fc Colafrancescol (|2003)). 
The range fluxes depends on the R e / P , the ratio of elec- 
trons to protons at relativistic energies. Using the same 
energy range and cluster temperature we find that sec- 
ondary electrons produce a flux of 1.9 x 10 -13 ergs cm~ 2 
s _1 . For our best fit photon index Fx — 2.21 and 90% 
confidence normalization range of 7.0-11.1 x 10 -5 pho- 
tons cm -2 s _1 keV -1 at 1 keV, we determine a flux of (6 
- 10) x 10" 13 ergs cm" 2 s _1 for the 0.13 - 100 keV range. 
This result has a lower limit 3 x greater than the results 
predicted for flux from secondary electrons. Although 
the 90% upper confidence limit is barely consistent with 
the primary electron model ( suggesting R e / P ~ 0.01), we 
point out that the results of IMiniati et al. I l)200UD is for 
shocks over the entire cluster, and our detection is only 
for a fraction of the cluster. That is the total non-thermal 
emission from cluster is probably greater than our detec- 
tion, further strengthening the argument that the flux is 
too high to come from secondary electrons. 
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